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Abstract 

We find supersaturation for the intrinsic gluon distribution of nuclei, i.e. the 
low x unintegrated nuclear gluon distribution peaks at intermediate trans- 
verse momenta kt = Q s and vanishes at zero kt- Taking into account the 
intrinsic transverse momenta of the gluons and the saturation scale p s of the 
produced gluons, we calculate the minijet cross section arising from gluon 
gluon scattering for RHIC energies. For central collisions at y/s = 200GeV 
the saturation scale p s ~ 1.4 GeV is larger due to intrinsic kt effects and 
increases with energy. Our theoretical results on charged particle multiplicity 

agree very well with the recent experimental data from RHIC. 
PACS number(s): 24.85.+p,25.75.+r,12.38.Bx 
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I. INTRODUCTION 



With the operation of RHIC |T|-0, heavy ion collisions have entered a new era, where 
semihard collisions become more important compared to soft collisions . Therefore a detailed 
knowledge of intrinsic parton distributions is important for a quantitative understanding 
of heavy quark production and/or cross sections at large transverse momenta. Due to 
the increase of the gluon distribution at small x, a sizeable cross section at mid rapidity 
comes from the region where the transverse momenta of partons become comparable to their 
longitudinal momenta. Therefore a k t factorization scheme is preferred, where the intrinsic 
transverse momenta of partons are included in the calculation of the minijet cross sections. 
Especially in nuclei, intrinsic transverse momenta of gluons can become quite large, since at 
small x the gluon clouds of the individual nucleons overlap and interact accumulating large 
transverse momenta. A careful inversion of the color dipole nuclear cross section allows to 
assess this effect. 

We differentiate between the gluon saturation scale Q s due to the interaction of gluons 
with one another inside the same nucleus before the collision and the saturation scale p s of 
produced gluons after the first collision between the nuclei. Both scales arise due to different 
mechanisms, but sometimes appear to be of the same magnitude. 

Inside the nucleus supersaturation occurs due to the color neutrality of the nucleons. 
One cannot have gluons at small x and small kt, because the nucleons from which they 
originate are color neutral. This finding is derived from a careful inversion of the multiple 
scattering of a test color dipole on the nucleus. It is confirmed by the evaluation of the 
BFKL equation in the nucleus [|J. We will give our derivation of the unintegrated gluon 
distribution in nuclei in section II. 

In section III we will use this distribution to calculate the mini jet cross section g + g — > 
g + g. Due to the intrinsic transverse momenta of the gluons the LO QCD cross section 
is enhanced at a given cut off scale. Since the intrinsic transverse momenta depend on the 
number of gluon clouds overlapping, we find that for more central collisions the enhancement 
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is larger. 

Section IV is devoted to the analysis of the gluon system which is produced after the first 
collision. These now on shell gluons are numerous at small nuclear impact parameters and 
crowd the impact parameter dependent area F(b) of intersection between the two nuclei. 
Concentrating on gluons in the central rapidity region one sees that their number ^ times 
their individual size vr/p^ cannot exceed the area F(b). This equation comes from an estimate 
of the number of merging gluons 2 — > 1 of produced gluons compared to the number of 
splitting gluons 1 — ► 2 |jU|. Since the number of gluons decreases severely as l/p 2 s and the 
saturation criterion F(b) xpj increases with p 2 s one obtains a self consistent determination of 
the saturation scale of produced partons p s . Intrinsic transverse momentum will increase 
the number of produced gluons at the same cut off. As a final result, these effects lead to 
a higher saturation scale p s ~ 1.5GeV than the scale found without intrinsic k t . In the 
framework of local parton hadron duality the so derived gluon distribution can be directly 
applied to estimate the number of produced charged hadrons. From e + e~ collisions || 
the conversion factor is known and of order unity. In this section, we will also present the 
prediction on charged particle multiplicity based on our saturation criterion. As first pointed 
in |9j , the centrality dependence of charged particle multiplicity can be used to test different 
saturation models. And very recently, there have been several studies on this quantity by 
using different saturation approaches |l0|,[ll| . We will find the saturation criterion we used 
in the following gives reasonable description on centrality dependence of charged particle 
reported by PHENIX collaboration [Q]. 

In section V we give a summary of our results and discuss the possibility to determine 
the intrinsic saturation scale Q s by measuring "monojets". 

II. UNINTEGRATED GLUON DISTRIBUTIONS IN COLD NUCLEI 

We know that the dipole-proton cross section depends on the unintegrated gluon distri- 
bution function f g (x,k 2 ). Assuming a perturbative 2 gluon exchange picture one has the 
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following equation for a color neutral dipole with transverse size r, virtuality Q 2 and energy 
s = [Q 2 + m 2 )/x. We assume a s fixed, fll2"f 

4.77- r rfih. 

**^(*,r) = -j J -^[l-e^]a s f g (x,k 2 ). (1) 

Here f g (x, k 2 ) is the unintegrated gluon distribution of the proton, which is related to the 
normal integrated gluon density as 

xG(x,Q 2 ) = f dk 2 f g (x,k 2 ). (2) 

The unintegrated gluon distribution function can be obtained from the color dipole- 
proton cross section &d-p m Eq. ([!]) by inversion. Since the cross section contains the color 
dipole factor [1 — e *], the inversion is more complicated than a Fourier Bessel transform. 
The energy of the dipole cross section under the integrand depends on the virtuality k\ as 
s = {k\ + m 2 )/x. 

a s f g (x,k 2 ) = ^kt^kt^f^JoM&^r) 

= / ^[-hrJxihr) - ^(Jofor) - J 2 {k t r))] x a d _ p (s,r), (3) 

l07T d J T 2 

where Jo, Ji, and J 2 are Bessel functions. 

We can derive the unintegrated gluon distribution function in nuclei by using the same 
method. The cross section of a test color dipole on the nucleus A depends on the nuclear 
thickness function T(b). For a hard sphere nucleus we would have Tib) = 2p ^J R\ — b 2 with 
Po = 0.16/m~ 3 . For the following calculations, however, we always use realistic thickness 
functions obtained from a nuclear density distribution with finite thickness. This is impor- 
tant for peripheral collisions. According to Glauber theory the total nuclear cross section 
has the following form: 

(4) 

a d - A {b, x, r) = 2(1 - e -hnb)^- P (s,r)y (5) 
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The differential cross section <Jd- A {b, x, r) with respect to the nuclear impact parameter 
b is related to the unintegrated gluon distribution f g / A (x,k t ,b) in the nucleus at impact 
parameter b. 

Air f cPk+ 

a d - A (b,x,r) = — J dPb-j^[l-e^]a a f g/A (x,kt,b), (6) 

which yields the normal gluon distribution in the nucleus after integration. 

rQ 2 

xG A (x,Q 2 ) = dk 2 d 2 bf g/A (x,k t ,b). (7) 

From Eqs. @ and we can derive the unintegrated gluon density of nucleus A as, 

a s f g/A (x,k t ,b) = ^J ^[-hrJ^hr) - ^f(J (k t r) - J 2 {k t r))} x 2(1 -e"^^-^). 

One sees that in the limit of small density the integrated density in the nucleus would 
be A times as big as in the proton (shadowing effects will give some modifications). To 
obtain the unintegrated gluon density of nuclei from the above equations, we must know 
the parameterization of the dipole cross section ad- p (r). Up to now, there are various 



parameterizations for this function |T3|-|T5|. As a rough estimate, we use the following 
simple parameterization, 

cr d - P (r) = C(0,s)r 2 . (9) 



We take the energy dependent coefficient C(0, s) from the parameterization of Ref. [15 



C(0,s) = a (s)/Rl(s), (10) 

where the parameterizations for cro(s) and Rq(s) can be found in ||15|| . To compensate other 
effects (e.g. shadowing |nj), we normalize our unintegrated gluon density to the integrated 
gluon density according to Eq. (|7|) for Q 2 ^> Q 2 . Finally, we get the following results for the 
double differential gluon distribution, (cf. Fig. 1) 

f g/A (x,k t ,b) = T{b) x x G A (x)-^-e~^, (11) 



where 



Q 2 = 2T(b)xC(0,s), (12) 

Note that the intrinsic momentum distribution depends on the closeness to the center of 
the nucleus. The smaller b the more gluon clouds in the nucleus overlap producing a larger 
gluon transverse momentum. At central 6 = the mean transverse momentum in the gluons 
increases with nuclear number A as (k 2 ) = 2Q 2 oc A^GeV 2 . The transverse momentum dis- 
tribution peaks at k 2 = Q 2 and then decreases towards small kt, vanishing at zero transverse 
momentum. We call this dynamical phenomena "supersaturation" . A similar behavior for 
the unintegrated gluon distribution function of nucleus was also derived by the evaluation 
of the BFKL equation in the nucleus [ffl. Supersaturation is quite different from saturation 



where growth of the gluon density slows down at small transverse momenta cf. ||17|| . The 
dynamical origin of supersaturation at small k t is the dipole-nucleus cross section for large 
dipoles, which is reduced due to multiple scattering. 

III. PRODUCED GLUONS AFTER ONE COLLISION IN HOT INTERACTION 

AREA 

Having the unintegrated gluon density of the nuclei, we calculate the production of 
gluons and quarks after one collision. To include the intrinsic k t effects, we adopt the 
phenomenological kt-kick model used in previous calculations of hard photon and hadron 
production |19|-21]. Following Ref. [20|, we write the cross section for hard collision as 



da(b, PT )_ r d 2 kTS{h . kT) d^A^ (13) 



dpr J dp' T 

where b is the impact parameter of the nucleus-nucleus collision. The differential cross section 
da jp, T ^ is identified as the standard LO QCD predictions for the 2 — > 2 processes without 
intrinsic k t effects, where the final produced parton would have the transverse momentum 
p' T . Due to the kick k t from the sum of the transverse momenta for the two incident partons, 
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the real final transverse momentum — p^ + fcr- The averaging function S(b; kx) reflects 
the intrinsic k t effects from the incident partons, which depend on the impact parameter b 
of the two nuclei, and can be obtained from the unintegrated gluon distribution functions 
calculated in the previous section. Since the intrinsic gluon distribution functions of the two 
nuclei depend separately on the locations b\ and b 2 of the event relative to each nuclear center, 
we first derive the differential averaging function depending on b\ and b 2 , from Eq. (|Tl~D , 



, ,11, , ,11,. ,:i 4t ,;i 
S{hMM = / ^^j^e-^j^e-^\q 1T + ^-k T ), (14) 

J 7T IT {Q sl ) Z [Qs2) 

where 61, b 2 , and b satisfy the relation bi — b 2 = b, and Q 2 sl1 Q 2 2 have the following forms, 

Qli = 2T A (h) x C(0, s), Q 2 s2 = 2T B (b 2 ) x C(0, s). (15) 

Since the rest of the integrand the differential mini jet cross section —j^- does not depend 
on b\, b 2 , we can integrate Eq. (0) over b\ and b 2 , and obtain a normalized averaging 
function only depending on b, 

o (h u ^_ /d 2 M 2 62TA(6i)T B (6 2 )g(6 1 ,6 2 ;fc r )^ 2 )(6 1 -6 2 -6) 

b ^* T) Jd*hd?b 2 T A fa)T B (b a )5W{h-b i -b) ' [ } 

The above averaging function S(b; k?) is normalized to unity after integration over k t . In 

Fig. 2, we plot this quantity as a function of k t for different impact parameters b. This figure 

shows that for more central collisions (smaller b), there is a stronger effect from intrinsic 

gluon transverse momenta k t . 

From the above formulas, we derive the number of gluons produced in the first collision 

of nuclei A + B at impact parameter b, applying an infrared cut off po on the produced 

gluons. 



where 



N A Ei{b) = T AB (b)a hard (b,p ), (17) 



T AB (b) = J d%T A (b - 6i)T B (&i), (18) 
a hard (b, Po )= I d p T d 2 k T S(b;k T )^4^. (19) 



We find that if we neglect the intrinsic k t effects (setting S = J( 2 )(for) in the above equa- 
tion), we return back to the results in QCD- improved parton model based on the collinear 
factorization approach [^]. However, for high energy nuclear collisions, we can not neglect 
these intrinsic k t effects, which enhance the cross section significantly at a given infrared cut 
off p . 



IV. SATURATION OF PRODUCED GLUONS AND INCLUSIVE CROSS 

SECTION 

For central collision (b = 0) of identical nuclei A+A with radius Ra, the interaction 
area of the two nuclei is identical to their transverse size ttR a . Assuming the radius of 
each produced gluon with transverse momentum po to be 1/po, we have saturation for an 
infrared cut off p$ = p s when the produced quanta fill the transverse interaction area irR\. 
All gluons with momenta smaller than the saturation momenta fuse into gluons with a total 
momentum larger than the saturation momentum. They are not realistic endproducts for 
the application of local parton hadron duality to the minijet cross section. So, for central 
collision, we have the following saturation condition, 

7T 7T IT R ^ 

N A a(Po = p s ,b = 0)- = T AA (0)a hard (p s )- = —±. (20) 

Pi pi (3 

The parameter (3 takes care of some modification of this simple geometrical picture which 
can be estimated from the evolution equation including gluon gluon fusion [g] or from || and 
is expected to be around 2 — 3. More realistic may be a simulation of classical Yang-Mills 
theory which gives (5 — 1.4 - 2.0 [p^,p^j . 



We extend the above saturation equation to any arbitrary impact parameter b of the two 
nuclei, 

7T 7T F(b) 

N A a(Po = Ps, b)— = T AA (b)a hard (p s (b); b)—^ = ——, (21) 
P s Psib) 2 p 

where F(b) is the interaction area for b > 0. From the geometry of the collision at arbitrary 
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impact parameter b and taking into account the realistic nuclear density profiles of two 
nuclei, we have the following form for the function F(b): 



where <7j„ is the inelastic cross section for proton-proton scattering. In our calculations, we 
use <r in = 35, 39, 45mb for y/s = 56, 130, 200GeV, respectively. In the limit of sharp-edged 
nuclear density, this area function will lead to the normal area function of two incident nuclei 
at impact parameter b 



Naturally, the saturation momentum becomes a function of the impact parameter p s = p s {b). 
Having determined the saturation scale p s (b) from Eq. ([H]), we plot in Fig. 3 the numerical 
results for p s as a function of b for three typical energies at BNL RHIC, y/s = 56, 130, and 
200GeV, where we set (3 = 1.8. From this figure, we can see that the saturation scale p s 
increases with collision energy. For central collision (b = 0), p s = 1.2GeV at y/s = 56GeV, 
increases to p s = 1.3GeV at y/s = 130GeV, and to p s = 1.35GeV at y/s = 200GeV. 

It is important to differentiate between the saturation scales Q s {b\),Q s (b-2) of intrinsic 
virtual gluons and the saturation scale p s (b) of produced gluons. The intrinsic saturation 
scale is a property of each individual cold nucleus and is defined at the local impact parameter 
in each nucleus. It parameterizes the intrinsic gluon distribution functions in both nuclei 
from which the averaging function S(b, fcr) is derived. The saturation scale for the produced 
gluons p s is a property of the interacting system of the two nuclei and depends on the nucleus 
nucleus impact parameter b. 

Our results show only a weak dependence of p s on b which is not so dramatic as found in 
ref. ||10|| . Furthermore, the saturation scale p s in our model is larger than lGeV for a wide 
range of b. We note that the intrinsic saturation scale Q s for the gluon distribution in each 
individual nucleus flT5|), however, is sensitive to the cylindrical distance b' of the center of 
the nucleus. In Fig. 3 we show both saturation scales for comparison. 



F(b) = J 



d 2 b lC l 2 b 2 {l - e - TAibl)a -){l - e- TA(bl) ° in )5 i2 \bx -b 2 - b), 



(22) 
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With the saturation scale p s (b) determined from Eq. (pl~D we can calculate the number 
of partons produced after the first collision of A + B, from which we can further estimate 
the charged particle multiplicity ||[/],|n]]. In Fig. 4 we plot the quantity / '(0.5N part .) 
as a function of N part .. The number of participants N part .(b) in an A+A collision can be 
calculated using 

N part . = 2 J d% 1 T A (b 1 )(l - e - T -^">). (23) 

The charged particle multiplicity is related to the parton number produced by the 
first collision = §^p- To take into account the difference of rapidity and the pseudo- 
rapidity, we include another factor 0.9 for ^ = 0.9^ > which can be calculated 
from above formulas. From Fig. 4, we see that our predictions on centrality dependence of 
charged particle multiplicity agree well with the PHENIX data as well as the PHOBOS data 
except few data points associated with very peripheral collisions, for which the saturation 
mechanism may not be valid any more. 

As a final remark, we note that in [|ItJ, quite different from our global saturation approach 
discussed above, a local saturation condition was introduced for every b\ and b% i n each 
nucleus. If the so determined local p s (bi,b 2 ) is getting too small, an additional cut off is 
implemented in ref. [|10[ such that p s (bi, 62) > 0.5 GeV. This local saturation approach leads 



to a stronger discrepancy for peripheral collisions than our model cf. Fig 4 and Fig. 4 of 
Ref. [[|. However, to finally settle this issue, we need more experimental observables and 
data as suggested in p3| . 



V. CONCLUSIONS 

In this paper, we have calculated the minijet production in A — A collisions at RHIC, 
including the enhancement from intrinsic transverse momenta of the incident gluons. We 
derived the unintegrated gluon distributions of cold nuclei and found "super saturation" ; the 
number density of gluons vanishes at zero transverse momentum and has a maximum at the 
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intrinsic saturation scale Q s . We differentiate between this gluon saturation scale Q s inside 
both nuclei before the collision and the saturation scale p s of produced gluons after the 
first collision between the nuclei. To obtain the saturation scale p s , we apply the saturation 
criterion in the central rapidity region, where the produced parton number ^ times their 
individual size Tr/p 2 s must not exceed the area F(b) of the two interacting nucleus. Intrinsic 
transverse momentum of the incident partons increases the number of produced gluons (i.e., 
the production cross section) at the same cut off, and then leads to a higher saturation scale 
p s pa 1.5GeV than the scale found without intrinsic k t . Finally, we find that our numerical 
results on the centrality dependence of the charged particle multiplicity agree well with the 
experimental PHENIX and PHOBOS data, which support the saturation scheme we used 
in this paper. 

We note as a final remark, that it is also very important to directly determine the 
intrinsic transverse momentum scale Q s for large nuclei. For this purpose we propose to 
measure "monojet" production at RHIC , in which the final state only contains one large 
Pt jet without any other large px balancing jets. By triggering on produced particles in the 
projectile and target fragmentation regions with a summed high transverse momentum in one 
azimuthal hemisphere, one may select minijets with transverse momentum in the opposite 
azimuthal hemisphere which come from 2 — >• 1 processes in the central region. In principle 
one can thereby test the intrinsic saturation scale Q s which is characteristic for isolated 
nuclei and also accessible in electron nucleus scattering. The same configuration of monojet 



production was also suggested in [pj] to search the jet quenching effects in high energy parton 
production. The difference between these two approaches and the experimental accessible 
will be discussed in the future studies. 
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Figure Captions 



FIG.l. The double differential gluon distribution f g /A(x,k t ,b) for Au(A = 197, Ra = 
6.37 fm) at x — 1CT 2 as a function of k t [GeV] at different cylindrical distances (b = 
0, 1/2Ra, 3/4Ra, Ra) from the central axis of the nucleus. The maxima of these dis- 
tributions are at k t = Q s « OMGeV, 0.90GeV, 0.76GeV, 0A6GeV for these four cases 
respectively. 

FIG. 2. The averaging function S(b; kr) vs kx for different nucleus-nucleus impact parame- 
ters b shows that intrinsic transverse momentum kicks k t with larger kr are more numerous 
for central collisions at b = than for peripheral collisions with large b. 
FIG. 3. The saturation scale p s in GeV is shown as a function of b for three different energies 
= 5QGeV, 130GeV, 200GeV. We also plot the intrinsic saturation scale Q s (x = 0.01) as 
a function of b' for the individual nuclei before collision, where b' is the cylindrical distance 
from the central axis of the nucleus. In contrast to p s , Q s is a more strongly varying function 
of b'. 

FIG. 4. The theoretical normalized charged particle multiplicity [dN ch /dr]]/ (0.5N part _) as a 
function of the number of participants N part . is shown for the theory with intrinsic gluon 
transverse momentum k t (full line) and without k t averaging effects (dashed line). The 
crosses mark the experimental PHENIX (+) and PHOBOS (x) data. 
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